Renormalization and destruction of I/7 2 tori in the standard nontwist map 
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Extending the work of del-Castillo-Negrete, Greene, and Morrison, Physica D 91, 1 (1996) and 
100, 311 (1997) on the standard nontwist map, the breakup of an invariant torus with winding num- 
ber equal to the inverse golden mean squared is studied. Improved numerical techniques provide 
the greater accuracy that is needed for this case. The new results are interpreted within the renor- 
malization group framework by constructing a renormalization operator on the space of commuting 
map pairs, and by studying the fixed points of the so constructed operator. 



In recent years, area-preserving maps that violate the twist condition locally in phase space have 
been the object of interest in several studies in physics and mathematics. These nontwist maps show 
up in a variety of physical models. An important problem from the physics point of view is the 
■ understanding of the breakup of invariant tori, which show remarkable resilience in the region where 
the twist condition is violated, called shearless tori. In terms of the physical system modelled, these 
tori represent transport barriers, and their breakup corresponds to the transition to global chaos. 
Mathematically, nontwist maps present a challenge since the standard proofs of celebrated theorems 
__, in the theory of area-preserving maps rely heavily on the twist condition. In this paper, we study the 
^"*^ ■ breakup of the shearless torus with winding number I/7 2 , where 7 is the golden mean. This torus 
serves as a test case for improved techniques we developed. At the point of breakup the shearless 



torus exhibits universal scaling behavior which leads to a renormalization group interpretation. 

I. INTRODUCTION 



In this paper we consider the standard nontwist map (SNM) M, as introduced in Ref. § 

£2 ' Xn+1 = Xn + a i 1 ~ Vn+l) 



Un+i = Vn ~ &sin(27ra„) , (1) 
where (x, y) € T x R, a € (0, 1), and b G (—00,00). The map M is area-preserving and violates the twist condition 

dx i+ i (xi,yi) 



t^O V(xi, yi ), (2) 



along a curve in phase space, which has been recently called the nonmonotone curve.^ Traditionally, most studies 
\ of area-preserving maps have dealt with the twist case, but in recent years more and more research has been focused 
on the nontwist case. 

Applications of nontwist maps occur in many fields, for example, the study of magnetic field lines in toroidal plasma 
devices (see e.g. Refs. |[ ||), in celestial mechanics, § fluid dynamicsjjj and atomic physics. |(| It has been shown |^| 
that nontwist regions appear generically in area-preserving maps that have a tripling bifurcation of an elliptic fixed 
point. In addition to these applications, the map is quite interesting from a mathematical standpoint because many 
important theorems in the theory of area-preserving maps assume the validity of the twist condition, e.g. the KAM 
theorem and the Poincarc-Birkhoff theorem. The SNM can serve as a model for the development of new proofs. Up 
to now, only a few mathematical results exist for nontwist maps (see ag. Ref. || |[ [ToJ pd[ ). 

We continue the work of del-Castillo-Negrete, Greene and Morrison, |l~2|, [l3| who studied the breakup of the shearless 
invariant torus with winding number I/7, where 7 = (1 + V5)/2 is the golden mean. We present the analysis of the 
breakup of the shearless invariant torus with winding number (in continued fraction representation) 

w= [0,2,1,1,...] = l/ 7 2 . (3) 

Because this winding number is a noble number (its continued fraction expansion ends with [1,1,1,...]), the behavior 
of the residues of the approximating periodic orbits is expected to be the same as in the I/7 case, i.e. we should find 
the same fixed point of the renormalization group operator with the same unstable eigenvalues that were found in 



Ref. 13L But, the form of the renormalization group operator, which is defined later in Sec. [V, is different from the 
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FIG. 1: Parameter space around the critical point (marked by *) of the uj — l/7 2 -shearless curve, showing the points for which 
shearless invariant tori exist. 



I/7 case. Also, the region of parameter space we study is different. Additionally, since the periods of approximating 
periodic orbits arc bigger than th ose fo r the 1 /j case, the present work serves as a test case for improved numerical 
techniques described later in Sec. Ill A . 

A different approach, which yields rough parameter values for the breakup of invariant tori, was used by Shinohara 
and Aizawa in Ref. who showed that a shearless invariant torus crosses the x-axis at two points, [fL5| xa = a/2 — 1/4 
and xb = a/2 + 1/4. For a given (a, 6), a point on the shearless torus, (a/2 + 1/4, 0), is iterated many times (we used 
10 6 ). If the y value stays below a threshold (we used \y\ < 0.52), it is assumed that the shearless curve exists and the 
point is plotted. Figure [l] depicts our duplication of their procedure. 

We see that the critical point for the I/7 2 shearless curve (indicated by *) lies on the boundary in Fig. [|. Thus, 
the boundary points of Fig. |l| represent the critical function for the SNM. This is a generalization of the definition of 
the critical function for the standard twist map (see e.g. Ref. |l6|), which has only one parameter e.g. k. The critical 
function in the twist case is then denned as k c (u>). Here, we have two parameters, but the shearless invariant torus of 
a given winding number ui exists only for parameter values belonging to a curve (a, b (a; a;)) in the parameter space. 
Thus, we can define the critical function by the critical points on each of those curves by (a c ,b(a c ]tu)). By finding 
the critical points for many other winding numbers (both nobles and non-nobles), we hope to find a more accurate 
critical function curve than the one shown in Fig. [l| 

In Sec. H we review some basic properties of the SNM. The detailed breakup of the shearless invariant torus with 
winding number I/7 2 is presented in Sec. [II, which also contains a discussion of the numerical procedures involved. 
In Sec. |v|, we interpret the results within the framework of the renormalization group. Section [v| a summary and 
some directions of future research. 



II. REVIEW OF RESIDUE CRITERION AND STANDARD NONTWIST MAP 



In this section, we give a brief review of some basic concepts of the the ory of area-preserving maps in the context 



of the SNM. For a more in-depth discussion the reader is referred to Ref. y and references therein 
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Periodic orbits and residue criterion 



Since the pioneering work of Greene, [pT| periodic orbits have proven to be very useful for studying the breakup of 
invariant tori in area-preserving maps. Below are some standard definitions. 

An orbit of an area-preserving map M is a sequence of points {(xi,yi)} such that M (xt,yi) = (aij+i, The 
winding number uj of an orbit is defined as uj = linii_ >00 (a; ! :/i) when the limit exists. Here the x-coordinate is "lifted" 
from T to R. A periodic orbit of period n is a sequence of n points {(xi, yi)}" =1 , such that M" (xi, j/j) = (xi + m, yi) 
for all i = 1, ... ,n, and m is an integer. Periodic orbits have rational winding numbers ui = m/n. An invariant torus 
is an orbit with irrational winding number that covers densely a one-dimensional set in phase space. Of particular 
importance are the invariant tori that wind around the x-domain because, in two-dimensional maps, they act as 
transport barriers. 
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The linear stability of a periodic orbit is determined by the value of its residue,\Yl\ R, which is defined as R := 
[2 — Tr(L)] /4. Here, L is the map M n linearized about the periodic orbit of interest and Tr denotes the trace. If 
< R < 1, the orbit is stable or elliptic; if R < or R > 1, it is unstable or hyperbolic; in the degenerate cases R = 
and R = 1, it is parabolic. 



Periodic orbits can be used to systematically approximate invariant tori. 1 17 The method is based on the observation 
that given a sequence of rational numbers {nii/rn} whose limit is lo, the sequence of periodic orbits with winding 
numbers {rrii/rii} approaches the invariant torus with winding number lo in phase space. It is important to find the 
"best" possible sequence, i.e. the sequence that converges to lo the fastest. The elements of the best possible sequence 
(see e.g. Rcf. |l|) are the convergents that arc obtained from successive truncations of the continued fraction expansion 

Of LO. 

The residue criterion |l7j] can be stated as follows: Consider an invariant torus with winding number lo. Let {mi/n{\ 
be the sequence of convergents approximating lo, and Ri the residues of their corresponding periodic orbits. 

1. If lim^oo \Ri\ =0, the invariant torus exists. 

2. If lirn^oo \Ri\ = oo, the invariant torus is destroyed. 

3. At the boundary in parameter space between those two limits, the invariant torus is at the threshold of destruc- 
tion and the residues either converge to a constant, non-zero value, or there are convergent subsequences. 

This criterion is based on the idea that the destruction of an invariant torus is caused by the de-stabilization of 
nearby periodic orbits. The residue criterion has been used successfully in many cases to predict with high precision 
the threshold for the destruction of invariant tori. Several theorems have been proved that lend mathematical support 
to the criterion. [[19, 20 



The numerical search for periodic orbits is difficult because, in principle, it is a two-dimensional root finding 
problem. However, the task is considerably simplified for reversible maps, |l) which are maps that can be factored 
as M = IioIq, where io,i are involution maps that satisfy l\ = I§ = 1. The sets of fixed points of the involution maps, 
To,! = {(x, y) |/o,i {x, y) = (x,y)}, are one-dimensional sets, called symmetry lines of the map. Once we kno w Tp,i 



the search for periodic orbits is reduced to a one-dimensional root finding problem, as explained below in Sec. Ill A 1 



B. Standard nontwist map 



The SNM is reversible. The symmetry lines Tq, composed of fixed points of Iq are si = {(x,y)\x = 0} and 
,S2 = {(x,y)\x = 1/2}. The symmetry lines T\, composed of fixed points of I\ are S3 = \[x,y)\x = a (l — y 2 ) /2} and 
s 4 = {(x,y)\x = a(l-y 2 ) /2 + 1/2}. 

A major difference between the standard nontwist map and twist maps is that there are two periodic orbits, if they 
exist, with the same winding number on each symmetry line. This can be seen easily in the integrable case. For 6 = 0, 
the m/n periodic orbits on the si symmetry line are located at 

{x, y) = (0, ±y/l - (m/n)/a) . (4) 

We will call the orbit with the bigger (smaller) y-coordinate the up (down) periodic orbit. 
The SNM is also invariant with respect to the transformation 

T(x,y) = (x+l/2,-y). (5) 

The coordinates of the up and down periodic orbits on the symmetry lines Si, denoted by (x U i,y U i) and (xdi,ydi) 
respectively, are related by this symmetry as follows: 

{x d2 ,yd2) = T((x ul ,y ul )) , (x u2 ,y U 2) = T((x dl ,y dl )) , 

{xdA,Vdi) = T((Xu3,yu3)) , {Xui,Vui) = T((xd3,yd3)) ■ (6) 

Therefore, it is actually enough to compute periodic orbits on si and S3, since the orbits along the other symmetry 
lines can be obtained from (p). 
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Periodic orbit collisions and bifurcation curves 



Periodic orbits in the SNM can undergo a particular form of bifurcation that occurs when the up and down 
periodic orbits of the same winding number meet on the symmetry line. These collisions were detected numerically 
in Rcfs. [22, 23 and [l^. Further studies of this bifurcation can be found in Refs. |^. [IT]. 

From (BR it follows that, for a given a. only periodic orbits with m/n < a exist at b = 0. As the value of b increases, 
the up and down orbits approach each other and at the bifurcation value, they collide and annihilate each other. For 
higher values of 6, both orbits no longer exist. Figure |^ illustrates the behavior of periodic orbits as we increase b 
from 6 = 0. Here the y-coordinates of the m/n = 3/8 periodic orbits on Si is shown as a function of b for the fixed 
value of a = 0.4. 



Up Orbit 
Down orbit 



0.3 0.4 

b 



FIG. 2: Plot of behavior of the up and down periodic orbit of winding number 3/8 for increasing 6-values at a = 0.4. The 
vertical axis shows the {/-coordinates of the orbits along si. 
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FIG. 3: Bifurcation curves for several convergents of 1/7 



Based on these numerical observations, the notion of a bifurcation curve in parameter space was defined in Ref. |l2j. 
The m/n-bifurcation curve b = <& m / n (a) is the set of (a, b) values for which the m/n up and down periodic orbits 
are at the point of collision. The main property of this curve is that for (a, b) values below b = $ m y„(a), the r/s 
periodic orbits, with r/s < m/n exist. Thus, m/n is the maximum winding number for parameter values along the 
m/n-bifurcation curve. 

The idea of approximating invariant tori with irrational winding numbers by periodic orbits is used to define the 
bifurcation curve for an invariant torus as follows :]l2[ The u) -bifurcation curve b = $ w (a) for an irrational ui is the 
set of (a, b) values such that b = $ w (a) = lim^oo <& m ./ n .(a), where <& m ./„.(a) is the TO.;/?ii-bifurcation curve and 
{mi/ni} are the convergents of w. For (a, b) points along the w-bifurcation curve the invariant torus with irrational 
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winding number ui is the curve of maximum winding number and is called shearless. Figure |3| depicts the bifurcation 
curves for several convergents of I/7 2 . This figure also makes it plausible that the limit in the above definition exists. 

III. BREAKUP OF TORUS WITH w = I/7 2 

In this section, we present the analysis of the breakup of the shearless invariant torus in the standard nontwist map 
with winding number uj = I/7 2 . Tables | and || list the convergents used for these calculations. For more details see 
Ref. |4j. 



w 


Fi/F i+ 2 


\i] 


Fi/F i+ 2 


[1] 


1/3 


[19] 


6765/17711 


[3] 


3/8 


[21] 


17711/46368 


[5] 


8/21 


[23] 


46368/121393 


[7] 


21/55 


[25] 


121393/317811 


[9] 


55/144 


[27] 


317811/832040 


[11] 


144/377 


[29] 


832040/2178309 


[13] 


377/987 


[31] 


2178309/5702887 


[15] 


987/2584 


[33] 


5702887/14930352 


[17] 


2584/6765 


[35] 


14930352/39088169 



TABLE I: Some of the convergents of to = [0, 2, 1,1,.. .] for which the periodic orbits still exist at criticality. 



w 


Fi 1 -fi+2 


[*1 


Fi/Fi+2 


[2] 


2/5 


[18] 


4181/10946 


[4] 


5/13 


[20] 


10946/28657 


[6] 


13/34 


[22] 


28657/75025 


[8] 


34/89 


[24] 


75025/196418 


[10] 


89/233 


[26] 


196418/514229 


[12] 


233/610 


[28] 


514229/1346269 


[14] 


610/1597 


[30] 


1346269/3524578 


[16] 


1597/4181 


[32] 


3524578/9227465 



TABLE II: Some of the convergents of u> = [0, 2, 1,1,.. .] for which the periodic orbits do not exist at criticality. 



A. Numerical methods 

The computational steps necessary to find the critical point and the residue behavior of the approximating periodic 
orbits are as follows: 

1. Find a good approximation to the l/7 2 -bifurcation curve in (a, 6)-space using the bifurcation curves for its 
convergents. 

2. Along this bifurcation curve, find the up and down periodic orbits on the symmetry line Si that approximate 
the invariant torus, and compute their residues. 

3. Locate the (a, b) point along the curve at which the residues exhibit critical behavior. 

4. Find the residues of the periodic orbits at criticality along the remaining symmetry lines. 

5. Find the eigenvalues of the unstable eigenmodes of the renormalization group operator. The details of how to 
do this depend crucially on the type of critical scaling behavior that is exhibited by the residues. 
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1. Searching for periodic orbits 

Periodic orbits on the symmetry lines can be be computed relatively easily for reversible maps using the following 
result :|§ If (x, y) G r 0) i then M n (x, y) = (x, y) if and only if M n / 2 (x, y) G T 0:1 (for n even) or M( n±1 )/ 2 (x, y) G r 1)0 
(for n odd). Thus, for example, periodic orbits with odd period n on the S\ symmetry line can be obtained by 
looking for points on s\ that are mapped to S3 or S4 after (n + l)/2 iterations. This can be implemented as a one- 
dimensional root finding problem by considering the zeros of the function F(y) = sin [2-7T (x — a (l — y 2 ) /2)] , where 
(x,y) := M(" +1 )/ 2 (0,y). The sine function is used to eliminate the difference between S3 and S4. Similar ideas can 
be applied to find other orbits. 



2. Finding m/n-bifurcation curves 

Recall that the bifurcation curve <& m / n (a) of a periodic orbit of winding number m/n was defined in Sec. |II C| to 
be the set of points (a, 6), at which the up and down periodic orbits of winding number m/n collide along the si 
symmetry line. Thus, at a given value of a, the function F(y) has two roots for b < $ m /„(a), no roots (locally) for 
b > 3> m / n (a) and a single root, which is also an extremum, for b = <& m / n (a). We thus search for the zero of the 
extremum of F(y) as b is varied. 

To find the whole (or large portions) of a bifurcation curve, we use the monotonic nature of the curve (see Fig. |^) 
as follows: Given a point (ai,6i) on the bifurcation curve i.e. b\ = $ m /„(ai), we increase a by a fixed amount to 
<22 = a\ + a st cp. We then start at the point (02, 61) and increase b until we reach 62 = <!>,„/„ (0,2). To make sure that 
we are finding the correct bifurcation curve, we start searching (a, 6)-spacc at (a, b) = (m/n, 0). Even then, the steps 
in a cannot be taken to be too large. Experience has shown that steps in a of 1 x 10~ 5 or 1 x 10~ 6 are safe. This 
method is unfortunately very slow because the part of the curve at small b values is very steep and the interesting 
(near critical) part of the curve is far away from the 6 = limit. 

We managed to drastically improve the speed of these calculations by using the following ideas: 

1. Numerical evidence strongly suggests that a bifurcation curve is smooth and monotonically increasing, although 
it is not proved mathematically. p5| So we use linear extrapolation from two previous points to find the new 
value of b around which to search for the bifurcation point. It was found that any higher order extrapolation 
did not improve the algorithm further. 

2. To find bifurcation curves for periodic orbits with very large periods (e.g. of the order of several million) the 
following procedure is used: Starting at the bifurcation curve of a smaller period, we increase b until the 
bifurcation curve of the higher period is reached. The advantage of this procedure is that we do not need to 
do the extremely time consuming calculations of the bifurcation curves for very high period orbits starting at 
6 = 0, but rather we can search for them near the region of interest. 



3. Finding 1 / '■y 2 -bifurcation curve 

Recall that the l/7 2 -bifurcation curve was defined as the limit of m.i/n.i bifurcation curves, where m,/nj are 
convergents of I/7 2 . It was numerically observed that close to criticality, this limit is approached in accordance with 
the following scaling relation: [p^t 

*[»+!] (o) = $ 1/f (fl) + B n (a) vl' V \ (7) 

where the $r n i(a) denotes the bifurcation curve of the periodic orbit with winding number [n] = F n /F n+ 2, v\ is a 
number to be determined later, and B n (a) is a period-twelve function, i.e. 5 n -|_i2 (a) = B n (a). 
If Eq. (^) holds, it follows that for fixed a 

$ 1/7 3 = lim , $ ["+l] $ ["+»]-^[n]^[~ + 13] (8) 
rwoc ($ [n+1] - $[„]) - ($ [n+ i 3] - $[„+12]j 

We obtained the l/7 2 -bifurcation curve using n = 19 in Eq. (g), i.e. using the bifurcation curves for [32], [31], [20] 
and [19] (see Tables | and ||). 

Now one can justify a posteriori the use of Eq. (@). Solving Eq. (§) with a = a c for v\ yields: 



n — >oo 



*[«+!] (Oc) 
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and 



B n (a c ) = ($[„+!] (a c ) - b c ) v x 



ij\2 



(10) 



where (a c , 6 C ) is the critical point for breakup of the shearless 1/7 invariant torus i.e. b c = 3>i/ 7 2 (a c ). We found that 
v 1 1 ^ 12 = 2.678. Some numerical evidence for the periodicity of B n (a c ) is given in Table III. 



n B n (a c ) 


n B n (a c ) 


15 -0.4865 

17 -0.7090 

18 0.5019 

19 -0.3901 


27 -0.4865 
29 -0.7078 
20 0.5028 
31 -0.3887 



TABLE III: Period-twelve behavior of the scaling function B n (a c ). 



B. Results 



In this subsection, we present the results of our computations. 



1. Residue behavior at criticality 

We computed bifurcation curves up to [32] = 3524578/9227465 and found the critical points along them, i.e. the 
parameter values along those curves for which the residues of approximating periodic orbits neither converge to zero 
nor diverge to infinity. Figure |I] shows the critical residue behavior of the up and down periodic orbits on the symmetry 
line S\ along several different bifurcation curves. For lower period bifurcation curves, the residues first show signs of 
a six-cycle (to be discussed later in greater details), but then converge to \Ri\ ~ 0.25. This is because the invariant 
torus we are studying is not quite shearless. Thus we see the same behavior of the residues as in the case of a twist 
map. As we proceed to higher period bifurcation curves, the behavior of the residues of the approximating periodic 
orbits found along the s\ symm etry line resembles more and more a six-cycle. A renormalization group interpretation 
of these results is given in Sec. |y| 

Finally, we found the critical point (a c , b c ) along the l/7 2 -bifurcation curve to be the following: 



= 0.425160543 



b r = 0.9244636470355. 



(11) 



At the critical parameter values (a c , b c ), the residues of the down periodic orbits on si, which are equal to the residues 
of the up periodic orbits on s 2 because of the symmetry of the map (see Eq. (|^)), converge to the six-cycle 26 
{Ci.Ca, 0,, Ci.C-B.Cfe}, where 



Ci = -0.609 ± 0.005, C 2 = -1.288 ± 0.002, 
C 3 = 2.593 ± 0.005, C 4 = 1.584 ± 0.008, 
C 5 = 2.336 ±0.006, C 6 = 2.593 ± 0.005. 



(12) 



The six-cycle can clearly be seen in Fig. |^ which shows the residues of the up and down periodic orbits at the critical 
point along the si symmetry line. The residues of the up periodic orbits on si (and of the down periodic orbits on 
s 2 ) converge to the six-cycle {Di,D 2 , D 3 , D4, D 5 , D 6 }, where 



£>i = 1.584 ±0.008, D 2 = -1.288 ±0.002, 
D 3 = -2.630 ± 0.006, D A = -0.609 ± 0.005, 
D 5 = 2.336 ±0.006, D 6 = -2.630 ± 0.006. 



(13) 



The residue convergence for other symmetry lines is shown in Table IV where we denote by R Ui and the residues 
of the up and down periodic orbit on the symmetry line s,;. Note that the six-cycle {Di} of R U1 and Rd 2 (respectively, 
the six-cycle {Ci} of R U2 and Rdt) is observed to be the same as that of R U3 and Rd 4 (respectively, R Ui and Rd 3 ) 




FIG. 4: Residue behavior of the up (top figure) and down (bottom figure) periodic orbits on si at the critical points on 
bifurcation curves of [20] (+), [24] (o)and [28] (*). 



9 




[n] [n] 

FIG. 5: Residue behavior of the up (left figure) and down (right figure) periodic orbits on si at the critical point on I/7 2 - 
bifurcation curve. 





Ru 1 — Rd,2 


Ru2 — Rd\ 


R113 — Rd^ 


R11.4 — 


[i] 


, [13] 


, [25] 


Dx 


Ci 


D 4 




[3] 


, [15] 


, [27] 


D 2 


c 2 


Dr, 


c 5 


[5] 


, [17] 


, [29] 


D 3 


c 3 


D 6 


c 6 


[7] 


, [19] 


, [31] 


Di 


c 4 


Di 


Ci 


[9] 


, [21] 


, [33] 


D 5 


c 5 


D 2 




[11] 


, [23] 


, [35] 


D 6 


c 6 


D 3 





TABLE IV: Period-six convergence pattern of the residues near criticality along the different symmetry lines. 



except it is shifted. The two six-cycles arc related because of the symmetry of the map as follows: D\ = C4, D 2 = C%, 
D4 = C\, D5 = C5, C3 = Cg, and D3 = D§, It was numerically observed that Cq w — Dq, and therefore C3 f» — D3. 
Using these relations we see that there are only five independent residues which we take to be Ci,C2,Cs,C4, and C5. 

We compared the values of the residues at three different points along the l/7 2 -bifurcation curve, one point be- 
low criticality, one at criticality, and one above criticality: (a_,6_) = (0.425160540,0.9244636195728), (a c ,b c ) = 
(0.425160543,0.9244636470355) and («+,&+) = (0.425160545,0.9244636653440), respectively. The numerical results 
for the Ci are listed in Table |V|. We see that each element of the six-cycle tends to zero for (a_ , 6_ ) , to infinity for 
(a+, &+), while it tends to the critical value at (a c , b c ). Figure [| clearly illustrates this behavior. 

A comparison with the results of Ref. [l2] shows that, within numerical accuracy, we found the same values for the 
residues of the six-cycle, but the sequence has shifted by two: C\ — H3, C2 = H4, C3 = H$, C4 = Hq, C5 = Hi and 
Ce — H 2 , where Hi denote the residues for the I/7 case found in Ref. [l2| A similar shift by two occurs for other 
symmetry lines. 



2. Spatial scaling at criticality 

As expected, the shearless curve exhibits scale invariance at criticality, which can be demonstrated explicitly by 
using symmetry line coordinates^^ {x,y) denned by x = x — a(l — y 2 )/2 and y = y — y s . In these coordinates, 
the S3 symmetry line becomes a straight line that intersects the shearless curve at the origin. We find that, in 
symmetry line coordinates, the shearless I/7 2 invariant torus at criticality remains invariant under a scale change 
(x, y) — > (a 12 cc, f3 12 y). This property is illustrated in Fig. 0. 
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M 




Z- 


Zc 


z+ 




Z- 


Zc 


z+ 


[01] 




0.565 


0.565 


0.565 




0.914 


0.914 


0.914 


[07] 




-0.702 


-0.702 


-0.702 




1.893 


1.893 


1.893 


[13] 




-0.601 


-0.601 


-0.601 




1.574 


1.574 


1.574 


[19] 




-0.611 


-0.611 


-0.612 




1.590 


1.591 


1.592 


[25] 




-0.610 


-0.612 


-0.614 




1.578 


1.591 


1.599 


[31] 




-0.566 


-0.612 


-0.644 




1.406 


1.581 


1.710 


[03] 


c 2 


-0.752 


-0.752 


-0.752 


c 5 


2.169 


2.169 


2.169 


[09] 




-1.328 


-1.328 


-1.328 




2.505 


2.505 


2.505 


[15] 




-1.286 


-1.286 


-1.286 




2.329 


2.329 


2.329 


[21] 




-1.289 


-1.290 


-1.291 




2.337 


2.340 


2.341 


[27] 




-1.273 


-1.289 


-1.300 




2.300 


2.338 


2.364 


[33] 




-1.161 


-1.249 


-1.276 




1.873 


2.288 


2.614 


[05] 


c 3 


3.450 


3.450 


3.450 


c 6 


3.450 


3.450 


3.450 


[11] 




2.534 


2.534 


2.534 




2.534 


2.534 


2.534 


[17] 




2.598 


2.598 


2.598 




2.598 


2.598 


2.598 


[23] 




2.588 


2.594 


2.598 




2.588 


2.594 


2.598 


[29] 




2.498 


2.588 


2.650 




2.498 


2.588 


2.650 



TABLE V: Numerical values of the residue six-cycle d at z_ = (a_,6_), z c = (<z c , 6 C ), and z+ = (a+,b+). 



As described in Ref. [13], we can find y s using 



lis 



lim 



V[2i+1] U[2i+n] ~ V[2i-1] 2/[2i+13] 



0.47253494777, 



(14) 



[V[2i+1] ~ V[2i-1]) - {V[2i+13] ~ V[2i+ll]) 

where denotes the y-coordinate of the periodic orbit [n] along the S3 symmetry line. To obtain the quoted value 
of y s we used i = 10. We then obtained a and j3 as follows: p7j 



and 



lim 



j3 = lim 



"[2i+l] 



x [2i+13] 



V[2i+1] 



2/[2i+13] 



1/12 



1/12 



1.61759 



1.65702, 



115) 



(16) 



where (xr n i,yr n ]) are symmetry line coordinates of the point of the periodic orbit [n] that is the closest to the origin. 
Within numerical accuracy, these values are the same as in Ref. [l3| 

Further numerical analysis shows that periodic orbits also exhibit scaling behavior locally near the S3 symmetry 
line. Figure [8] shows points of the periodic orbit [n] = [21] (in symmetry line coordinates) and points of the periodic 
orbit [33] with the x and y coordinates rescaled by a 12 and /3 12 respectively. The result suggests that periodic orbits 
remain invariant under a simultaneous spatial rescaling and shifting of the winding number by twelve from [n] to 
[n + 12]. 



C. Numerical accuracy 



We conclude this section with comments about the numerical accuracy of the results. 

1. Points on the m/n-bifurcation curves were found with an accuracy ranging between 10 -12 and 10 -15 , where the 
larger value corresponds to larger periods. We obtain this measure of accuracy from the condition 

F(y) = 0, F'(y)=0 and F"(y) ? 0, (17) 

as explained previously. The numbers quoted above are the values of F(y) obtained at the numerically found 
minima in y. 
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FIG. 6: Residue Convergence for Ci through Ce at Zi (+), Z2 (o) and 23 (x) (see Table ^). 

2. Periodic orbits along the different symmetry lines around the critical point where found with an accuracy ranging 
between 10 -15 and 10 -17 . Here, the criterion is the difference between the winding number m/n of the periodic 
orbit of interest, and the winding number of the orbit that results when starting at the numerically found 
location of the periodic orbit on the respective symmetry line, and then iterating the map. 

3. A criterion for the accuracy of the scaled bifurcation curve $i/ 7 2 is harder to find, since the location of the 
actual curve is unknown. An upper bound on the error, though, should be the distance between the $i/ 7 2 
and $[37], since the latter definitely lies on the other side of $i/ 7 2. This error was found to be approximately 
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FIG. 7: Invariance under rescaling of shearless I/7 2 torus at criticality. 
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FIG. 8: Invariance of periodic orbits under simultaneous rescaling and shift of winding numbers by twelve. Here we show the 
periodic orbits [21] (x) and [33] (o) after rescaling of x and y coordinates by a 12 and (3 12 respectively. 



2 X 1CT 13 . 

4. A criterion for the accuracy of the critical point in parameter space, (a c , b c ) is even harder to define, since we 
cannot actually find the residues of existing orbits of all periods, which is required to check if the six-cycle at 
that point continues ad infinitum. We believe that the value for a c is accurate up to 1 x 1CP 9 . 

5. The uncertainties for the critical residues, quoted above in Eq. (|l2"|), were computed from the variation in 
numerical values of residues at a c for the three or four highest period orbits found (e.g. for C±, using residues 



13 



of [13], [19], [25] and [31]). 

6. If we evaluate the residues of the up and down periodic orbits along the Si and S3 symmetry lines, then the 



residues on the other two symmetry lines can be constructed using symmetry arguments (see Sec. II A). But, 
as a check of the numerical procedures, wc independently evaluated the residues on all the four symmetry lines 
and confirmed the symmetry arguments. 

IV. RENORMALIZATION GROUP INTERPRETATION 

In this section, we interpret the above results within the renormalization group framework. The analysis follows 
Refs. |lj closely, since, as expected, the residue behavior exhibits a six-cycle at criticality. But, because of the different 
winding number (i.e. different from I/7) the renormalization group operator will have a different form. 

Renormalization ideas have been used fruitfully in area-preserving maps and Hamiltonian flows. (See e.g. Refs. Eq, 



|29| , 30, BlL |32|, [33[ |3j| and references therein.) In contrast to mathematical KAM theory, which proves the existence 
of dense sets of invariant tori in regions of phase space, the renormalization group approach addresses the problem of 
the destruction of an invariant torus with a specific winding number under strong perturbation. 



The following renormalization approach (see e.g. Refs. 13, 3^) is based on the residue criterion (Sec. II A). To study 



the breakup of an invariant torus of winding number u>, we loosely represent the map M as: 

M= {R 1 ,R 2 ,R 3 ,...), (18) 

where the {Ri} are the residues of the periodic orbits with winding numbers {rrii/rii}, the convergents of w. For 
example, the integrable map will be represented by (0, 0,0,.. .) because all the orbits are parabolic in that case. The 
key idea is to construct an operator 1Z that explores the infinite tail of (|l^) by mapping a map given by ( |l8| ) to another 

map, 1Z(M), represented by 1Z (AI) = (Ai, R2, ■ ■ ■ , Ri, ■ ■ ■ j , where Ri = R4+1. This operation can be interpreted as a 

time renormalization since periodic orbits with large periods are transformed into periodic orbits with smaller periods, 
which amounts to a rescaling of time. 

The residue criterion can now be rephrased in this framework: 

1. If lim n ^oo TZ n (M) = (0, 0,0,.. .), the invariant torus exists. 

2. If linin^oo lZ n (M) = (±00, ±00, . . .), the invariant torus is destroyed. 

3. If limn^oo 7Z n (M) is a map for which the residues have finite, nonzero values, i.e. a map that is invariant under 
the action of lZ m (a fixed point of lZ m ) for some m > 0, the invariant torus is at the threshold of destruction. 
Possible scenarios are the convergence of the residues to a fixed value or to a convergent subsequence. 

There are two kinds of fixed points: simple fixed points and critical fixed points. In the case of area-preserving maps, 
we come to the following interpretation. A simple fixed point is an integrable map (all the residues are zero), and 
its basin of attraction contains all the maps for which the invariant torus exists. A critical fixed point is a map for 
which the invariant torus under consideration is at criticality. All the maps in its basin of attraction exhibit the same 
universal behavior at the critical breakup. 

A. Renormalization group operator 

Following the discussion in Refs. ^ and [l3|, we use pairs of commuting maps because they provide a simple way to 
define the renormalization operators for invariant tori. 

A pair of commuting maps is an ordered pair of maps, (U, T), such that UT = TU. An orbit of a point (x,y) 
generated by (U, T) is the set of points {U m T n {x, y)}, where m and n are integers. A periodic orbit of period m/n is 
an orbit for which U m T n (xi,yi) = (xi,yi). 

For the breakup of the invariant torus with winding number u> = I/7 2 , we define the renormalization group operator 
by 





(19) 

As for the case of 1/7,|13| this operator contains both time and space renormalization as follows: 
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The space renormalization is represented by the operator B, which rcscalcs the (x,y) coordinates, i.e. (x, y) — > 
B(x, y) where 

At the critical fixed point studied in this paper, we see that r = a and s = (3 given by Eqs. (|ll|)-(|l(f). 

The time renormalization is, again, accomplished by the specific combination of the commuting maps. If (x, y) is 
a periodic orbit of (U,T) with winding number Fi/Fi + 2, then B(x,y) is a periodic orbit of (U,T) = 1Z(U,T) with 
winding number Fi^i/ F i+ i, as can be verified as follows: 

U Fi -^T Fi +^B{x,y) = B{UT)- F ^ (UT 2 ) F ' +1 (x,y) 

= BU- F ^ +F ^ T- F ^ +2F ^(x,y) 
= BU F > T F+2 {x,y) 
= B(x,y). 

By induction, an orbit with winding number Fi/Fi + 2 under (U,T) is transformed into an orbit of lZ n (U,T) with 
winding number Fi- n /F i+2 - n . 



B. Simple periodic orbit of 7Z 

We can find the intcgrable period-two orbit (U±,T±) of the renormalization operator (|l9| ) by requiring that 
7t(U±,T±) — (U T ,T T ). This two-cycle is given by the following pairs of maps: 



U± 



x 



x~i 2 t y 2 h 2 
y 



T, 



x 



x + 1 ± y 2 
V 



(21) 



where the rescaling of the coordinates is given by 



B 



-7 
±7 



(22) 



Using the definition U±T±(x,y) = (x,y) of the periodic orbits of period m/n, we get the rotation number as a 
function of y: 



w±{y) 



l±y 2 



l 



-i 2 Ty 2 h 2 7 



r 



2\ - 1 



1± 1 



y 



Thus we see that the map (C/-, T_) is locally equivalent, under a change of coordinates, to the SNM with parameters 
(a,6) = (l/ 7 2 ,0). 



C. Critical periodic orbit of 1Z 



The next step is to analyze the critical periodic orbit of 1Z. Consider the nontwist map 

C = (Ci, — , C 2 , — , C 3 , — , C4, — , C 5 , — , C 6 , — , Ci, — , C 2 , . . .) , (23) 

where the Ci are the elements of the six-cycle computed earlier, and the "-" denote the periodic orbits that do not 
exist (see Table |J). By construction, this map is a period-12 orbit of the renormalization group operator (a fixed 
point of 1Z 12 ), i.e. 

K 12 C=C. (24) 
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In Sec. IIIB, we found that the residues of the periodic orbits approximating the l/7 2 -shearless curve in the standard 
nontwist map exhibit convergence to the six-cycle {Cj}. Assuming that we can fine-tune the results for (a c , b c ), we 
expect that lim^^ TZ n M (a c , b c ) = C. 

If we are studying the breakup of the l/7 2 -shearless curve for parameter values along the bifurcation curve for 
one of the low-period convergents, then we start near the stable manifold of the critical periodic orbit of 1Z. But, 
under the action of 7Z, we are pushed along an unstable direction. Thus, we see parts of the six-cycle of residues (see 
Fig. |^), but the limiting residue behavior is observed to be lim^oo \Ri\ w 0.25, which is characteristic for the critical 
fixed point of twist maps (see e.g. Rcf. |33| ). In renormalization group language, this means that part of the unstable 
manifold of the critical nontwist fixed point (maps for which (a, b) is below (a c , 6 C )) is in the basin of attraction of the 
critical twist fixed point. 



D. Eigenvalues 



As shown in Ref. [T^, it is possible to use our numerical data to draw further conclusions about the renormalization 
group operator 7Z, in particular to compute its unstable eigenvalues. The main difficulty in computing these eigenvalues 
is that the space of maps is infinite-dimensional whereas the (a, b) parameter space has only two dimensions. But 
the fact that we can find an isolated point (a c ,b c ) in parameter space at which the map is at criticality means that 
the dimension of the unstable manifold is two. The map M at (a c ,6 c ) is the intersection point of the two-parameter 
family of maps with the stable manifold of the fixed point, i.e. values a c and b c describe the location of the critical 
fixed point of 1Z 12 in its unstable manifold. 

Below, we first compute the eigenvalues characterizing the approach to (a c , b c ) in the (a, b) parameter space using 
the numerical results from above. As shown in Ref. [l3| based on the type of numerical data obtained, the two 
eigenvalues can be found by 



^ =M (T7Tr) (25) 

n->oo y $[ n ] [a c j -b c J 

and 

V2 = Um ( <M2«+12]-°A (26) 
n^oo y a c [ 2n ] - a c J 

The last step is to relate the values i>i to the unstable eigenvalues Si of the renormalization group operator 7Z. The 
key idea is to study the behavior of the residues of the periodic orbits approximating lu at the (a, b) values used in 
the computation of For details see Ref. p"3[ 

Denoting the unstable eigenvalues of 1Z by S\ and 62, we conclude that Si = (1/vi) 1 ^ 12 . We find numerical values 

of 



61 as 2.678, S 2 « 1.583. (27) 

Comparing this to the values found in Ref. [l^ shows that within (assumed) numerical uncertainty these values are 
the same as those for I/7, as predicted. The a c values used to determine 82 were a c [ 2 6] and <J C [14]; which explains the 
larger discrepancy. Work is under way to improve this result. 



V. CONCLUSION 



We have shown through numerical simulations that the critical residue values at the breakup of the 1 /7 2 -shearlcss 
curve in the standard nontwist map coincide with those of the l/7-shearless curve. In addition, the critical scaling 
parameters and the unstable eigenvalues of the renormalization group operator were found to be the same for both 
cases. The main differences are the location of the respective critical point in parameter space and the detailed form 
of the renormalization group operator in terms of commuting maps pairs. 

Future work includes the search for the breakup values of more winding numbers to map out the details of the 
critical function depicted in Fig. [l]. In addition, new fixed points of the renormalization group operator might be 
obtained by studying the breakup of shearless curves with non-noble winding numbers. 
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